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We present a search for long-lived doubly-charged Higgs bosons with signatures of high 

ionization energy loss and muon-like penetration. We use 292 pb~^ of data collected in pp collisions 
at ^/s = 1.96 TeV by the CDF II detector at the Fermilab Tevatron. Observing no evidence of 
long-lived doubly-charged particle production, we exclude and bosons with masses below 
133 GeV/c^ and 109 GeV/c^ , respectively. In the degenerate case we exclude mass below 

146 GeV/c^ . All limits are quoted at the 95% conhdence level. 
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The electroweak gauge symmetry of the standard model (SM) is broken by the hypothetical Higgs mechanism, 
thereby imparting masses to the W and Z bosons, the mediators of the weak force. A number of models [1-4] extend 
the SM Higgs sector to include additional symmetries. For instance, the left-right symmetric model [2] postulates a 
right-handed version of the weak interaction, whose gauge symmetry is spontaneously broken at a high mass scale, 
leading to the parity-violating SM. This model is supported by recent data on neutrino oscillations [5], and explains 
small neutrino masses [6]. The model generally requires a Higgs triplet containing a doubly-charged Higgs boson 
which could be light in the minimal supersymmetric left-right model [3,4]. Discovery of the boson 

would not only shed light on the Higgs mechanism, but also provide evidence for new symmetries beyond the SM. 
Grand unified theories containing Higgs triplets and their relevance for neutrino masses and mixing are reviewed in [7], 
while “Little Higgs” models that ameliorate the heirarchy and fine-tuning problems of the SM are reviewed in [8]. 

The dominant production mode at the Tevatron is pp 1 * jZ + X ^ + X, whose cross section at 

tree level is specified by the quantum numbers and the mass (m//±±) of the boson. The partial width in the 

leptonic decay modes is given by Fo' = hfi,m where hui are phenomenological couplings. In a previous 
Letter [9], we published the most stringent mass limits from direct searches in the ee, ep and pp decay channels for 
0.5 > hill > 10“®. In this Letter, we discuss the case where the boson lifetime (r) is long (ct > 3m, corresponding 
to hii' < 10“®), resulting in the boson decaying outside the CDF detector [10]. A supersymmetric left-right 
model [4] has predicted a light boson with B — L = 0, where B and L represent baryon number and lepton 

number respectively, resulting in hui = 0 and a long lifetime [9]. The LEP experiments have set limits on a long-lived 
boson [11,12], with the best limit coming from the DELPHI experiment [12], excluding mH±± < 99.6 GeV/c^ 
(99.3 GeV/c^ ) at the 95% confidence level (C.L.) for bosons with couplings to left- (right-)handed leptons. Our 
search for pair-production of long-lived, doubly-charged particles is based on the signatures of increased ionization 
energy loss and muon-like penetration of shielding (due to their large mass). We set the most stringent mass 

limits in the context of the left-right symmetric model. 

This analysis uses 292 ± 18 pb“^ of data collected by the CDF H detector [13] in pp collisions at ^/s = 1.96 TeV at 
the Tevatron. The detector consists of a cylindrical magnetic spectrometer with silicon and drift chamber trackers, 
surrounded by a time-of-flight system, pre-shower detectors, electromagnetic (EM) and hadronic calorimeters, and 
muon detectors. The central drift chamber (COT) [14], central calorimeter [15] and the muon detectors [16], covering 
the region \'q\ < 1 [17], are used in this analysis. The COT and calorimeter provide ionization information in addition 
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to tracking and identification of penetrating particles. 

We use an inclusive muon trigger requiring a COT track with transverse momentum px > 18 GeV/c [17], and a 
matching track segment in the central muon chambers. In the offline analysis, we search for H pair-production 
by requiring two COT tracks, each with px > 20 GeV/c, beam impact parameter < 2 mm and at least 30 (out of a 
maximum of 96) sense wire hits. At least one of the tracks is required to have a matching muon chamber segment. We 
also require their isolation /o ,4 < 0.1, where / 0.4 is the ratio of the total calorimeter Ax [17] around the track within 
a cone of radius R = yj (Ap)^ -p {A(j))‘^ = 0.4 to the track px [17]. Energy deposited by the particle is excluded from 
the calculation of Iqa- Finally, we tag and reject cosmic ray tracks using an algorithm based on COT hit-timing [18], 
whose efficiency is measured to be 100 lg'g% for collider muons and leaves negligible cosmic ray contamination. 

We use Z —>■ p/x events that were triggered by one of the muons to measure trigger and offline identification 
efficiencies of the other muon. The track selection efficiency is (93.6 ± 0.2)%, and the efficiency for one of the two 
bosons to satisfy the muon trigger and matching-segment requirements is (96.8 ± 0.7)%. The effect of increased 
multiple-scattering of doubly-charged particles is investigated by comparing the segment matching efficiency for muons 
from Z boson decays with that for lower-px muons from T decays. The small («0.5%) difference, when scaled as 
to the large px of tracks, predicts a negligible (« 0 . 2 %) correction. About 3% of particles are expected 
to be sufficiently slow {(} < 0.4) to have a reduced efficiency due to delayed hits, for a net efficiency loss of 0.4%. 
A correction is applied to the track selection efficiency for bosons passing near a calorimeter tower edge and 

depositing a large ionization energy signal in an adjacent tower. This effect, caused by the resolution of the track 
extrapolation, leads to the boson candidate failing the isolation requirement. This geometrical correction results 
in an overall track selection efficiency of (89 ± 4)%. 

The charge collected by each COT wire is proportional to the ionization deposited by the particle per unit length 
(dA/dx), and is encoded in the width of the digital pulse generated by the front-end electronics [14]. Offline corrections 
are applied for the electronics response, track polar angle, COT high voltage, drift distance, drift direction with respect 
to track direction, gas pressure, attenuation along the sense wire, radial location of the sense wire, and time. The 
mean number of hits on our selected tracks is 85. The mean {w) of the lower 80% of the corrected widths of all 
recorded hits of a track is used as a measure of its ionization energy loss. The use of the truncated mean reduces the 
sensitivity to Landau fluctuations. 

The most probable dA/dx for a minimum-ionizing particle corresponds to w « 15 ns, as seen from the cosmic-ray 



muon distribution in Fig. 1. For the search we require w > 35 ns. The w distribution of the latter is modelled by 
quadrupling the w measurements of cosmic ray muons, as given by the (charge)^-dependence of ionization energy loss 
in the Bethe-Bloch equation. We use low-momentum protons from secondary interactions to measure the efficiency 
of the dB/dx cut on tracks, which are expected to have similar or greater dE/dx than said protons (see 

Fig. 1). We obtain a control sample enriched in highly-ionizing protons by selecting low-momentum positively-charged 
secondary [19] tracks. The pion contribution is statistically removed by subtracting the w distribution of negatively- 
charged secondary tracks. Using the resulting w distribution of protons, we measure the w selection efficiency to be 
> 99.5%. 



FIG. 1. The distribution of the COT dE/dx variable w for positively-charged secondary [19] tracks in the momentum range 
of 300 — 350 MeV/c (solid), for high—px cosmic ray muons (dashed), and the expectation for tracks (dotted). The latter 
is modelled by quadrupling the w measurements of cosmic ray muons. The arrow indicates the signal selection region. 


We perform two simultaneous searches with “loose” and “tight” selections for highly-ionizing particles. The loose 
selection, based on the COT dE/dx measurement only, yields the maximum acceptance, while the tight selection 
also requires large EM and hadronic calorimeter signals for confirmation of a potential signal. We make the a priori 
decision to use the results from the “loose” search to quote an upper limit on the signal cross section, and the “tight” 
search results to quote a statistically significant observation of signal. The most probable ionization energy signal 
deposited by muons in the EM and hadronic calorimeters (referred to as Eem and Ehad respectively) is 0.3 GeV and 
1.7 GeV respectively, for normal incidence. For the tight search we require Eem > 0.6 GeV and Uhad > 4 GeV. 
The efficiency of the calorimeter ionization requirements is (81.1 ± 0.1)%, measured by quadrupling Eem and Uhad of 
a pure cosmic ray sample to model the energy deposition. 

We calculate the geometric and kinematic acceptance for a pair of bosons using the PYTHIA [20] generator 
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TABLE I. Summary of fake rate measurements. The e, p and r fake rates and the “muon fake rates” for jets are quoted as 
upper limits at the 68% C.L., since no events in the respective control samples pass the selection cuts. 


source 

“track” 

loose search 

“muon” 

“track” 

tight search 

“muon” 

jet (xlO'"^) 

O Q-|-O.U 
0.^_2.9 

< 0.05 

U.ZO_Q Q5 


< 0.05 

e (xlO-®) 

< 4 

< 0.00009 

< 0.05 


< 0.00002 

p(xlO-®) 

< 7 

< 7 

< 0.02 


< 0.02 

r (xl0“®) 

< 2 

< 0.002 

< 2 


< 0.002 


and a GEANT [21]-based detector simulation. The acceptance increases from 38.4% at mH±± = 90 GeV/c^ to 
46.8% at mH±± = 160 GeV/c^ , with the dominant relative systematic uncertainty of 1% due to parton distribution 
functions (PDFs) [22]. Systematic uncertainties due to momentum scale and resolution are negligible. 

Backgrounds arise from (1) jets fragmenting into high-px tracks, (2) Z ee, (3) Z /rp, and (4) Z ^ tt where at 
least one r decays hadronically. The backgrounds are a result of muon misidentification and dE/dx mismeasurement, 
which can arise from overlapping particles. Each background is estimated by multiplying the number of misidentifiable 
events by the product of the appropriate misidentification probabilities (fake rates). Fake rates are measured with 
and without the requirement of a matching muon chamber segment. We refer to these as the “muon fake rate” 
and “track fake rate”, respectively. A fake rate is defined as the probability that a track (or muon) passing certain 
loose identification cuts also satisfies the analysis cuts. For jets, electrons and r’s, the muon fake rate is obtained by 
multiplying the track fake rate by the estimated probability of mis-matching a muon chamber segment to the track. 

The track fake rate and muon fake rate for jets are measured from jet-triggered data and muon-triggered data, 
respectively. The variation of the fake rates with px and jet proximity is taken as a measure of systematic uncertainty. 
The number of misidentifiable jet events is given by the number of muon-triggered data events containing a loosely- 
selected muon and another loosely-selected track. Fake rates for electrons and hadronically decaying r’s are estimated 
from the GEANT-based detector simulation. These fake rate measurements are limited by Monte Carlo statistics, as 
no Monte Carlo events pass the selection cuts. The number of misidentifiable Z ^ ee events is obtained from the 
Z ee data sample, corrected for electron efficiencies and normalized to the luminosity of the muon-triggered signal 
sample. The number of Z —>■ tt misidentifiable events is obtained from the number of Z pp events observed in the 
data, assuming p — r universality, and correcting for muon efficiencies. Finally, fake rates for muons are measured 
from a pure sample of cosmic rays, which are again statistically limited as no events pass the selection cuts. The 
number of misidentifiable events is given by the number of Z fj,fx data events selected with the loose cuts. Table I 
summarizes the fake rate measurements, and Table II summarizes the resulting background estimates. 
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TABLE II. Summary of the estimated number of background events (quoted as 68% C.L. upper limits) and the observed 
number of events in the data. 


background 

loose search 

tight search 

jet 

< 3 X 

< 3 X lO-*^ 

Z ^ ee 

< 1 X 10““ 

< 2 X 10““ 

Z —> fijj, 

< 4 X 10“’^ 

< 4 X 10““ 

Z —> TT 

< 8 X 10“® 

< 8 X 10“® 

data 

0 

0 


No candidate events are found in the data. The null result is used to set upper limits on the number 

of signal events (3.2 at the 95% C.L.) and the pair production cross section using a Bayesian [23] approach, 
with a flat prior for the signal cross section and Gaussian priors for the uncertainties on acceptance, background 
and integrated luminosity (6%) [24]. The 95% C.L. upper limit on the cross section (which varies from 39.7 fb at 
mH±± = 90 GeV/c^ to 32.6 fb at mH±± = 160 GeV/c^ , see Fig. 2) is converted into an mass limit by 
comparing to the theoretical pp ^ "f* /Z X ^ + X cross section at next-to-leading order [25] using the 

CTEQ6 [22] set of PDFs. We include uncertainties in the theoretical cross sections due to PDFs (5%) [22] and higher- 
order QCD corrections (7.5%) [25] in the extraction of the mass limit. The theoretical cross sections are computed 
separately for 77^^ and 77^^ bosons that couple to left- and right-handed particles respectively. When only one of 
these states is accessible, we exclude the long-lived 77^^ boson below a mass of 133 GeV/c^ and the long-lived 77^^ 
boson below a mass of 109 GeV/c^ , both at the 95 % G.L. When the two states are degenerate in mass, we exclude 
mH±± < 146 GeV/c^ at the 95 % C.L. 



FIG. 2. The comparison of the experimental cross section upper limit with the theoretical next-to-leading order cross sec¬ 
tion [25] for pair production of 77** bosons. The theoretical cross sections are computed separately for bosons with left-handed 
(77)[ ) and right-handed (77)^*) couplings, and summed for the case that their masses are degenerate. 


In conclusion, we have searched for long-lived doubly-charged particles using their signatures of high ionization and 
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muon-like penetration. No evidence is found for pair-production of such particles, and we set the individual lower 
limits of 133 GeV/c^ and 109 GeV/c^ , respectively, on the masses of long-lived and bosons. The mass 
limit for the case of degenerate and bosons is 146 GeV/c^ . 

We thank M. Miihlleitner and M. Spira for calculating the next-to-leading order production cross section. 

We thank the Fermilab staff and the technical staffs of the participating institutions for their vital contributions. 
This work was supported by the U.S. Department of Energy and National Science Foundation; the Italian Istituto 
Nazionale di Fisica Nucleare; the Ministry of Education, Culture, Sports, Science and Technology of Japan; the 
Natural Sciences and Engineering Research Council of Canada; the National Science Council of the Republic of 
China; the Swiss National Science Foundation; the A.P. Sloan Foundation; the Bundesministerium fuer Bildung und 
Forschung, Germany; the Korean Science and Engineering Foundation and the Korean Research Foundation; the 
Particle Physics and Astronomy Research Council and the Royal Society, UK; the Russian Foundation for Basic 
Research; the Comision Interministerial de Ciencia y Tecnologia, Spain; and in part by the European Community’s 
Human Potential Programme under contract HPRN-CT-2002-00292, Probe for New Physics. 


[1] T. P. Cheng and L.-F. Li, Phys. Rev. D22, 2860 (1980). 

[2] R. N. Mohapatra and J. C. Pati, Phys. Rev. Dll, 566 (1975).; G. Senjanovic and R. N. Mohapatra, Phys. Rev. D12, 1502 
(1975); R. N. Mohapatra and G. Senjanovic, Phys. Rev. D23, 165 (1981). 

[3] C. S. Aulakh, A. Melfo, and G. Senjanovic, Phys. Rev. D57, 4174 (1998); Z. Chacko and R. N. Mohapatra, Phys. Rev. 
D58, 015003 (1998). 

[4] C. S. Aulakh, K. Benakli, and G. Senjanovic, Phys. Rev. Lett. 79 , 2188 (1997). 

[5] Super-Kamiokande Collaboration, Y. Ashie et ai, Phys. Rev. Lett. 93, 101801 (2004) and references therein. 

[6] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[7] B. Bajc, G. Senjanovic and F. Vissani, Phys. Rev. D 70 , 093002 (2004) and references therein; S. F. King, Rept. Prog. Phys. 
67 , 107 (2004) and references therein. 

[8] T. Han, H. E. Logan, B. McElrath and L.-T. Wang, Phys. Rev. D 67 , 095004 (2003). 


12 




[9] CDF Collaboration, D. Acosta et ai, Phys. Rev. Lett. 93, 221802 (2004). 

[10] The range 10“® > hui > 10“® corresponds to the boson decaying within the detector, and requires other triggering 
and tracking methods. 

[11] OPAL Collaboration, G. Abbiendi et al, Phys. Lett. B 526, 221 (2002). 

[12] DELPHI Collaboration, J. Abdallah et al., Phys. Lett. B 552, 127 (2003). 

[13] D. Acosta et al., Phys. Rev. D 71, 032001 (2005). 

[14] T. Affolder et al., Nucl. Instrum. Meth. Phys. Res. A 526, 249 (2004). 

[15] CDF Collaboration, F. Abe et ai, Nucl. Instrum. Meth. Phys. Res. A 271, 387 (1988). 

[16] G. Ascoli et al., Nucl. Instrum. Meth. Phys. Res. A 268, 33 (1988). 

[17] CDF uses a cylindrical coordinate system in which (j> is the azimuthal angle, r is the radius from the nominal beamline, and 
+z points in the direction of the proton beam and is zero at the center of the detector. The pseudorapidity rj = —ln[tan(0/2)], 
where 9 is the polar angle with respect to the 2 axis. Calorimeter energy (track momentum) measured transverse to the 
beam is denoted as Et (pt). Track px is evaluated from its curvature assuming a singly-charged particle, and is half of 
the true px for a doubly-charged particle. 

[18] A. V. Kotwal, H. K. Gerberich and C. Hays, Nucl. Instrum. Meth. Phys. Res. A 506, 110 (2003). 

[19] “Secondary” tracks are those produced in secondary interactions with detector material, and are selected by requiring 
beam impact parameter > 1 cm. 

[20] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994), version 6.127. 

[21] R. Brun and F. Carminati, CERN Program Library Long Writeup, W5013, 1993 (unpublished), version 3.15. 

[22] J. Pumplin et al., JHEP 0207, 012 (2002). 

[23] K. Hagiwara et al. (Particle Data Group), Phys. Rev. D 66, 010001 (2002); 1. Bertram et al., Fermilab-TM-2104, April 
2000 (unpublished). The numerical error caused by the upper cutoff in the numerical integration of the posterior likelihood 
is 0.005%. 

[24] S. Klimenko, J. Konigsberg, and T. M. Liss, Fermilab-FN-0741 (unpublished). 

[25] M. Miihlleitner and M. Spira, Phys. Rev. D68, 117701 (2003). The cross sections have theoretical uncertainties of (5-10)% 


13 



due to QCD corrections. 


14 



